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1 Introduction 



Non-perturbative phenomena of QCD such as confinement can be explored using 
Monte Carlo simulations of lattice gauge theory. The potential V{R) between two 
static quarks at separation R in quenched QCD is one of the simplest manifesta- 
tions of confinement and has been studied intensively. At large R the potential 
rises linearly as predicted by the hadronic string model. One can also measure 
the spatial distribution of the colour fields around such static quarks in order 
to get a detailed picture of the confining fiux tube. In Ref. [jll, which contains 
references to earlier work, this was done for the ground state and the first ex- 
cited state of the two-quark potential, having the symmetries of the Aig and Eu 
representations, respectively, of the lattice symmetry group D^h- Transverse and 
longitudinal profiles of chromoelectric and -magnetic fields were compared with 
vibrating string and dual QCD models for the fiux tube, with the latter model 
reproducing quite well the shape of the energy profile measured on a lattice. In- 
stead of SU(3), the gauge group used was SU(2), which is more manageable with 
present-day computer resources and is expected to have very similar features of 
confinement. This is reflected in the fact that the flux-tube models considered 
do not distinguish between SU(2) and SU(3). 

The method used to study the colour flelds on a lattice is to measure the cor- 
relation of a plaquette □ = |Tr(l — Un) with the Wilson loop W{R,T) that 
represents the static quark and antiquark at separation R. When the plaquette 
is located at t = T/2 in the /i, u plane, the following expression isolates, in the 
limit T oo, the contribution of the colour fleld at position r: 



In the naive continuum limit these contributions are related to the mean squared 
fluctuation of the Minkowski colour flelds by 



'^Bl{v) with J, k cyclic and /]^(r) ^ --E^{v). (2) 



When the interquark separation axis is chosen as the 1-axis the squares of the 
longitudinal and transverse electric and magnetic flelds can be identifled as 
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These can then be combined naively to give the action density 
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and the energy density 
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of the gluon field. 

Since in this work we use a plaquette to probe the colour fiux, the spatial size 
of the probe will decrease as the lattice spacing a ^ 0. To define a continuum 
limit of the colour fiux distributions, one would have to use a probe of a fixed 
physical size as a — >■ 0. In this work we wish to compare fiux distributions at 
different lattice spacing. One special tool that is available, when a plaquette 
is used to probe the colour fiux with the Wilson gauge action, is that exact 
identities can be derived for the integrals over all space of the fiux distributions. 
These sum rules Q relate spatial sums of the colour fields measured using 
Eq. m to the static potential V{R) via generalised /3-functions, which show how 
the bare couplings of the theory vary with the generalised lattice spacings in 
four directions. One can think of these sum rules as providing the appropriate 
anomalous dimension for the colour fiux sums. This normalises the colour fiux 
and provides a guide for comparing colour fiux distributions measured at different 
a-values. The full set of sum rules allow these generalised /5-functions to be 
determined at just one /3-value and references therein]. Here we investigate 
this further by comparing estimates at two different /5-values. This can also help 
to clear up some inconsistencies between the /5-function determination from the 
sum rules at one /5-value [|] and other methods p, |^ . 

In Ref. [0 the simulations were carried out at /9 = 2.4 with a 16^ x 32 lattice. Here 
results from similar computations at = 2.5 with a 24^ x 32 lattice are reported. 
Most of the simulation and analysis techniques, such as use of a variational basis 
with different fuzzing levels, are the same and can be found from Ref. [|l|. In 
addition to more accurate measurements of flux-tube proflles, we also present 
estimates of /5-functions at both /3 = 2.4 and 2.5. 

In Sect. 2 static quark-quark potentials V{R) are extracted and in Sect. 3 the 
corresponding flux tube proflles calculated. In Sect. 4, V^(-R) and the proflles are 
related by sum rules and various estimates of the /^-function are made. Some 
conclusions are made in Sect. 5. 

2 Static potentials 

We construct lattice operators to create and destroy states with two static quarks 
at separation R joined by gluonic paths which represent the colour flux. The 
techniques we use to make efficient operators with a large overlap with the ground 
state are described in detail in Ref. In order to improve the signal we applied 
"fuzzing" , where each spatial link is replaced by a weighted sum of itself and its 
surrounding spatial staples, before the correlations were measured. To investigate 
gluonic excitations and minimize their contribution to the ground state signal we 
need a variational basis, which was obtained by performing the measurements on 
lattices with different levels of fuzzing. 
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At /? = 2.4, 2.5 fuzzing levels 40, 16, and 40, 13, 2, respectively, formed the 
variational basis in the case of paths with Aig symmetry. A three-state basis 
may be expected to give a reasonable signal for the A'lg excitation by reducing 
the contamination from higher excitations with this symmetry. For paths with Eu 
symmetry the fuzzing levels 16, 13 were used for the two /9's with two different 
transverse extents of the paths forming a variational basis. These transverse 
extents were one and two lattice spacings for small longitudinal lengths and one 
and four lattice spacings for larger i?'s. In order to set the scale, the lattice steps 
a(2.4) and a(2.5) were determined by fitting the two-body parameterisation 

" e ■ 



ViR) 



+ bsR + Vo 
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to measured potentials at interquark separations R = 2, 3, 4, 6, 8, R = 2, 3, 4, 6, 12 



at /3 = 2.4, 2.5 respectively. Here 



is the latticized form of the Coulomb 



potential 1/r due to one-gluon exchange. The above range of R was chosen to 
correspond to similar physical distance ranges. The usual method for estimating 
the lattice spacing is to equate the dimensionless value of 65 from the fit of 
Eq. ^ to an experimental value. This is equivalent to utilizing V{R) in the limit 
i? — s> 00. However, our experimental knowledge of the two-quark potential comes 
from heavy mesons with r.m.s. radii around 1 fm. An alternative method due to 
Sommer compares the force from V{R) to experimental values at a distance 
range more appropriate to these mesons, i.e. r ^ 0.5 fm. In practice the equation 
(RoYF^Rq) = c is used to find Rq, where F(r) is the force between two static 
quarks at separation r. Various nonrelativistic continuum potential models give 
aRo ~ 0.49 fm for c = 1.65 and aRo ^ 0.66 fm for c = 2.44 The resulting 
scales and ratios of lattice spacings p = a(2.4)/a(2.5) are shown in Table |l|. The 
values of p from 65 and Sommer's scheme are seen to agree. 



p 


c= 1.65 


c = 2.44 


= 440 MeV 


2.4 


0.1098(5) 


0.1183(5) 


0.1190(5) 


2.5 


0.0778(4) 


0.0839(4) 


0.0846(4) 


p = a(2.4)/a(2.5) 


1.412(14) 


1.410(13) 


1.406(11) 



Table 1: Lattice spacing a determined using different methods 



The (3 = 2.5 potentials above have the accurate interpolations 

V{R)a^, = 0.555 + 0.0343i? - 0.280/i? (7) 
V{R)e^-V{R)a,, = 3.8/i? - 12.6//?^ + 24.6/i?^ - 18.8/i?^ (8) 
V{R)a[^ - V{R)au_, = 5M/R - 23.6/i?2 + 49.2/ R^ - 38.8/i?^ (9) 

The relations are valid for 2 < i? < 12, with no physical interpretation intended 
for the V{R)eu and V{R)a' expressions. The values per d.o.f. are 0.09, 1.35 
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and 0.22 for the Aig, and A[g states respectively. The corresponding potential 
fits for /? = 2.4 are given in Ref. [|l|]. 

Our variational basis is constructed from fuzzed link operators that represent 
creation or annihilation of two quarks at separation R, with the colour field in a 
specific state of lattice symmetry. These operators can be expanded in terms of 
the eigenstates of the transfer matrix 

1^) = co|l^o) + ci|\/i) + ... (10) 

with the measured correlation of a generalised Wilson loop given by 

W{R,T) = {Ro\Rt) = cle-""'^ (l + h{T/2y + . . .) (11) 

where 

h(t) = %-(y^-yo)t . (12) 

Co 

To minimize excited state contamination we need h <^ 1. 

As plaquettes in the middle of the generalised Wilson loop in the time direction 
are used to probe the colour fields, the relevant estimate of contamination is 
taken at t = T/2. The measured correlation is {Ro\Df\Rrp) instead of Eq. |ll[ 
This produces off-diagonal terms such as (Vi|nt|Vo) which increase the coefficient 
of excited state contribution from h to 2h for T = 2, t = 1. From the generalised 
Wilson loop ratios at each i?-value, we define V{T) = — ln[W{T)/W{T — 1)] since 
its rate of approach to a plateau as T — »• oo enables us to estimate the excited 
state contamination to the ground state. We calculate h from 

m = r/2)| ^ ^^V(T - 1) - V(T) = A ^'^ - ^> - 7 °°> . (13) 

Here the T — > oo extrapolated potential is defined as 

V{T ^ oo) ^ V{T) - X ^ e"(^-^»). 

1 — A 

In practice A was calculated from potentials at T = 1. Table ^ shows the excited 
state contamination for states of the two-body potential at /? = 2.5. It is seen 
that the Aig states are reasonably pure with the and A[g states containing 
increasing amounts of contamination. Note that the method used to extract the 
higher excited state contamination h to the first excited Aig state is unreliable 
since ground state contributions can dominate in principle. 



3 Colour field distributions 

Before embarking on the extraction of /3-functions using sum rules as presented 
in the next section, it is of interest to check further the quality of the input 
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Table 2: Excited state contamination at f3 = 2.5 as measured by h 



R 


state 


t = l 


2 


3 


2 




0.012 





005 


- 


3 




0.027 





012 


0.009 


4 


Aig 


0.036 





016 


0.008 


6 




0.073 





034 


0.016 


12 


Al9 


0.162 





075 


0.006 


2 


TP 
J^u 


0.2o2 





152 


O.lOo 


3 


Eu 


0.253 





131 


0.105 


4 


Eu 


0.250 





126 


0.076 


6 


Eu 


0.255 





123 


0.045 


12 


Eu 


0.427 





263 


0.189 


2 




0.201 





080 


0.076 


3 


A' 


0.228 





098 


0.070 


4 


A' 


0.261 





116 


0.122 


6 


A' 


0.333 





155 


0.157 


12 


A' 

^Ig 


0.518 





241 


0.204 



data for the sum rules by testing its scaling properties between [3 = 2.4 and 2.5. 
Essentially the input data is of two distinct forms - the two quark potential V 
and the spatial sums of the colour fields. In the past the scaling properties of V 
have been confirmed many times and will not be repeated here. However, scaling 
of the different colour field combinations is less clear, since the colour fields are 
measured using a plaquette, whose physical size changes with (3. This is relevant 
because only observables with the same physical size at different values of the 
coupling have a continuum limit. In the case of the three-dimensional sums over 
the colour fields, the lattice sum rules provide the appropriate normalisation as 
a — * 0. Even in this case, the scaling behaviour is only known after the divergent 
self-energies are subtracted. However, other observables, like the two-dimensional 
sums over transverse planes or the transverse profiles of the flux-tube, do not have 
a well defined scaling behaviour, but it is still of interest to explore how similar 
the profiles are at the two values of coupling used. 

From the relations in Eq. |^ and sum-rules to be presented in Eqs. p!^-p!6|, we 
can see that the measured action sums must be multiplied by the anomalous 
dimension h/ (3 to get the physical value, while the energy sums have a correction 
/ which goes to one in the continuum limit. Here 6, / are generalised /3-functions 
to be discussed later. As will be seen in sect. ^ below, the differences in the 
values of h and / between (3 = 2.4 and 2.5 are sufficiently small to neglect them 
in the following plots; e.g. 6(2.4)/6(2.5) = 0.97(7). Thus the normalisation of the 



5 



overall three dimensional sum over colour flux can be treated as almost constant 
in our j3 range. In turn this implies that we should compare more differential 
distributions using this scale. 

3.1 Spatial sums 

In Figs. the three-dimensional spatial sums of the action 5* and longitudinal 
and transverse energies El, Et involved in the sum rules are plotted as a function 
of R and T for the flux-tube ground state Aig and the two excited states Eu and 
A'^g for (3 = 2.5 and the scaled data at /3 = 2.4. The basic data {S, B) are 
dimensionless and require the factor to give energy and action densities in 
GeV/fm^. Since Figs. show the volume integrals of the basic data with (3 = 
2.5, for a scaling comparison, the volume integrals of the basic data for (3 = 2.4 
are multiplied by 2.4/(2.5p), where p is the ratio of lattice spacings given in Table 
|l|. The resulting two sets of data [P = 2.4(scaled), 2.5] should not be expected to 
lie on top of each other, since they have self-energies and self-actions that diverge 
as g'^/a from one-gluon exchange in leading order perturbation theory. However, 
the two sets of data should be parallel to each other, since the self-energies are 
independent of R. In the following extraction of the beta-function, it is the extent 
to which the slopes of these lines are non-zero that is relevant. 

In Figs. for clarity, only the data for one or two T values are drawn. The 
data for higher T have larger errors, but in most cases they are consistent with 
the data shown - indicating that the plateau in T is achieved. In Fig. |l|, for 
the action, the Aig and E^ states show scaling and with distinct non-zero slopes. 
This is best seen for the Aig state and deteriorates progressively in going to the 
Eu and A[g states. Scaling is questionable for the latter. Even so, the action 
and its i?-dependence are comparable in all three cases. A very approximate 
estimate of the difference ASo in the (3 = 2.5 and (3 = 2.4(scaled) self-actions is 
given by the vertical difference between the two sets of data. The result given by 
linear flts for the Aig ground state, A^o ~ 0.12(2), is not inconsistent with the 
Eu case, whereas the A'l^ excitation has a positive self-action difference only for 
the smallest R's. This curious feature of the A[g data is perhaps not surprising, 
since already Table § shows that this state has considerable contamination from 
other states. However, in spite of this, we thought it useful to include such data 
in this paper to illustrate where the data need to be improved. 

In Fig. H, for the longitudinal energy, the data are an order of magnitude smaller 
than for the action. Here the dependence on j3 is less clear and the presence of a 
non-zero slope much less distinct for the E^ and A[g cases. The difference A£^o 
in the f3 = 2.5 and f3 = 2.4(scaled) self-energies is now best taken visually at the 
lowest i?'s, giving AEq ^ 0.020(5) for the Aig state. The data for the gluonic 
excitations are consistent with this estimate. For Aig the curves cross at -R = 6 
and at higher i?'s the f3 = 2.4 curve lies higher. 



6 



In Fig. ^, for the transverse energy, only the Aig state shows a slope that is 
slightly non-zero, with both the and A[g cases having slopes consistent with 
zero. However, there is a distinct self-energy effect with AEq ~ 0.04(1) for the 
Aig state, again consistent with the gluonic excitations. 




Figure 1: The scaling of action (S) summed over the spatial lattice for paths with 
Aig, Eu and A'^g symmetries. 

Since it is the presence of non-zero slopes that is relevant for extracting the beta- 
function. Figs. and |^ already indicate that problems will arise when attempting 
to utilize the El^t data. 

3.2 Transverse sums 

The self-energy differences in the above spatial sums can be seen much more 
clearly in the transverse sums shown in Figs. ^|-^, where the longitudinal depen- 
dence of the sum over is presented for paths with Aig, E^ and A'^g symmetries 
and interquark separation i? = 8, 12 at /5 = 2.4, 2.5 respectively. In these fig- 
ures, i?L = corresponds to the mid-point of the interquark separation, while 
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Figure 2: The scaling of longitudinal energy {El) summed over the spatial lattice 
for paths with Aig, Eu and A'^g symmetries. 

Rl = R/2 corresponds to the position of the static quark sources. The R values 
were chosen to correspond to approximately the same physical distance at these 
two couplings, namely 0.946(4) and 1.007(5) fm. The plotted data is taken at 
T = 3, where we have a good signal to noise ratio. Unfortunately this means 
that the excited state contamination is relevant at T = 1, 2, where h is largest as 
can be seen from table ^ Furthermore, we are also using the largest -R's, where 
the excited state contamination becomes quite significant especially for the E^ 
and A[g cases. 

These transverse sums do not, strictly speaking, have a continuum limit. How- 
ever, in string models the transverse sums near the centre of long strings should 
be independent of R. So that, to the extent that string models are applicable 
and that R is sufficiently large, scaling would be expected. This is the assump- 
tion made in presenting the data in earlier works 0. In the figures, the basic 
(3 = 2.5 data is compared with the basic (3 = 2.4 data which has been multiplied 
by 2.4/(2.5p^). The longitudinal energy is plotted at half-integer lattice spacings 
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Figure 3: The scaling of transverse energy (Et) summed over the spatial lattice 
for paths with Aig, Eu and A'^g symmetries. 



(i.e. averaging over neighbouring values of £t, instead of £l, Bt) to get a bet- 
ter determination of the self-energy peak. In this case these peaks are expected 
to diverge as g"^ /a^ in physical units. 

Fig. ^ shows for the Aig state that, within error bars, all three field combinations 
(action. El, Et) scale well near the centre of the flux tube {Rl = 0). However, 
near the quarks [Rl ~ 6 in this case) the (3 = 2.4 and 2.5 curves differ consid- 
erably. This is mainly the effect of the self-energies being seen. This data near 
Rl = 6 shows some interesting features: 

a) Except for the f3 = 2.4 action, all the data exhibit a distinct peak near R^ = 6 
and the ASq, AEq extracted from Figs. 1-3 are qualitatively consistent with 
the corresponding values estimated from the areas under these peaks. We can 
estimate the self-energy peak height by subtracting the value at i?L = from the 
value at Rl ~ 6. For the action sums at 2.4 we cannot see any self-energy peak, 
while for El the ratio of self-energy peak heights at (3 = 2.4 and 2.5 is 0.59. For 
Et the peak height ratio is somewhat lower at 0.49; this will be discussed below. 



9 



o P=2.4 Action 
p=2.4 E^xlO 
|3=2.4 E^xlO 
|3=2.5 Action 
p=2.5 E^xlO 
p=2.5 E^xlO 




Figure 4: The dependence on longitudinal position (Rl) of the sum over the 
transverse plane of the colour flux contributions corresponding to the action, 
longitudinal (El) and transverse energy (Et) sum rules, eqs. Here Rl 

is measured from the mid-point for separation i? = 8, 12 at /3 = 2.4, 2.5 re- 
spectively. The data are in units of a(2.5) for the symmetric ground state (Ai^ 
representation) at T = 3. 

In comparison, the perturbative expectation (2.5a(2.5)2)/(2.4a(2.4)2) = 0.52(1) 
is in the middle of these values. This suggests that these peaks are dominated 
by self-energy effects. 

b) The transverse energy Et is completely dominated by the self-energy with 
the latter being, at least an order of magnitude larger than the non-self-energy 
terms which are expected to be essentially independent of Rl between the two 
quarks. This immediately explains the small slope of the Aig curve in Fig. ^ 
It also shows that any volume integral of the self-energy contribution cannot be 
accurately evaluated on the present lattice, since - in the f3 = 2.5 case - the 
Rl = 6 contribution is an order of magnitude larger than those from Rl = 5, 7 
i.e. the whole volume integral is given by contributions (each with ^ 10% errors) 
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Figure 5: As in Fig. ^ but for the first gluonic excitation {Eu representation). 
For each data set one error bar is shown; others are similar. 

from only three values of R^. Therefore, when the volume integrals from different 
values of R are subtracted, it will be hard to get a meaningful signal for the in- 
dependence of interest. 

In principle, the self-energy contribution can stretch out to the midpoint between 
the two quarks. According to our data, for the largest interquark separations as 
shown here, any such contribution seems negligible. This is seen by looking at 
the data near R^ ~ 10 — 12, which should be dominated by any self-energy tail. 
However, for the smaller i?'s the transverse sums at the midpoint can well have 
significant self-energy contributions, which is seen in later in section |4.3| when 
this data is used to determine /5-functions. 

c) The trend is that, for the action, the peak is about a 5% effect compared with 
the plateau contribution from Ri = 0, . . . , 6, for the peak is about 50% of 
the plateau contribution and for Et the peak completely dominates. Therefore, 
it should be expected that any predictions that depend on cancelling the self- 
energies are most reliable for the action and least for Et- 
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Figure 6: As in Fig. ^but for the second gluonic excitation {A[g representation). 
For each data set one error bar is shown; others are similar. 

In Fig. 1^ the corresponding transverse sums are shown for the state. The 
most notable features are: 

a) The action and El both scale within error bars for Rl ^ 0. 

b) As expected, the (3 = 2.4 data shows, for ^ 0, an enhancement of Et over 
its Aig counterpart, whereas the action and El are comparable to the Aig data. 
However, this enhancement in Et is not seen for j3 = 2.5, but this could be due 
to the relatively large errors for this case. 

c) Again the action for f3 = 2.4 does not exhibit a distinct peak near Rl = 6. 
Instead it simply shows a monotonic increase as Rl goes from to 6. The ratio 
of self-energy peak heights is 0.61, 0.59 for the action. El respectively - the latter 
being the same as in the Aig case. For E^ the ratio is 0.46, lower than for El 
by a similar amount that was observed for Aig. As the self-energies should be 
isotropic, the difference in the peak height ratios at the two /3's is worth exploring 
further. 
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This difference is caused by the different discretisation of longitudinal and trans- 
verse plaquettes. The electric field is dominant, so Et is mainly composed of 
plaquettes lying in the transverse plane, whereas El is mainly a planar sum of 
plaquettes with a perpendicular orientation. A simple way to investigate their 
diverging behaviour is to consider a scalar field (j) = e"^''"' with h ~ 10/fm 
and the integrals over it analogous to our sums. The (normalised) integral 
over the transverse two-dimensional plane, where the source lies, is found to 
change faster when the lattice spacing is varied than the integral over a three- 
dimensional "slice" of width a, analogous to E^. For our a(2.4) and a(2.5) the 
ratio [Plane(2.4)/Plane(2.5)]/[Slice(2.4)/Slice(2.5)] ^ 0.8 (also for = e"^!'^!'/^^), 
which agrees with the corresponding ratio observed for the Et and E^ for paths 
with Aig and E^ symmetries. The lower peak height ratios for Ex and the high 
Y^Et peaks observed in Figs. are also due to this effect. 

The diverging of the transverse sums over colour sources is quite consistent with 
the expectation from leading order perturbation theory. El diverges slower than 
the perturbative expectation, as expected, because of the transverse extent of 
the dominant electric field. On the other hand, for Et we would expect g'^/a? 
behavior, while the observed peak height ratios diverge slightly faster than this. 

d) Compared with 5* and El-, Et has a self-energy that is comparable to or larger 
than the plateau contribution from i?^ = 0, . . . , 6. 

In Fig. ^ the data for the A'^g state is shown with the following features: 

a) The action and Et are approximately scaling within the rather large error 
bars. However, for El (not shown) it is not possible to make this claim since the 
error bars are too large. 

b) The Et data, unlike that in the E^ case, now exhibits some enhancement - 
compared with the Aig state - for both f3 =2.4 and 2.5. For example, at f3 = 2.4 
and Rl = 0, ^T(^ig) ~ 0.01(1), ^t(^«) ~ 0.04(1) and ^T(^'ig) ~ 0.09(5). 
Even so, the plateau terms are still, at most, only comparable to the self-energies. 
Therefore, as for the Aig and Eu cases, those predictions that require a delicate 
cancellation of the self-energies are possibly not reliable. 

3.3 Transverse profiles 

In Figs. |l|-§| the scaling properties of the 3-d volume and, assuming the flux 
tubes have string-like features, of the 2-d transverse integrals are demonstrated 
- with some combinations of the colour fields being more successful than others 
in satisfying this property. It is, therefore, of interest to proceed finally to the 
"scaling" properties of the individual flux tube profiles. The transverse depen- 
dence of the action (S), longitudinal and transverse energies (El^t) ^ measured 
at the midpoint {Rl = in our convention) with the separations = 8, 12 at 
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j3 = 2.4, 2.5 respectively - is presented in Figs. |^-^. The correlations shown were 
measured at T = 3 for the Wilson loop. As we are again using a small T and 
large R, excited state contamination is significant in the Eu and A'lg cases. Here 
the (3 = 2.4 data is compared with the (3 = 2.5 data by multiplying the former 
by 2.4/(2.5p4). 

Aside from any intrinsic non-scaling arising from the different scale of the pla- 
quette used to probe the flux distributions, we should also be aware that effects 
can arise from the discretisation versus Rt of the distribution and, possibly, from 
self-energy effects. We have found earlier that the self-energy effects are negligible 
at the midpoint {Rl = 0) in the integrated distributions. We here assume that 
this applies to the differential distributions so this contribution can be neglected. 
The effect of the discretisation in R^ is that a sharply peaked distribution will be 
suppressed at coarser lattice spacing. There is some sign of this latter effect in 
our data: The smaller plaquette at j3 = 2.5 should increase the observed height 
of peaks (such as the centre of the flux-tube Rt = in the Aig case), whereas 
at the smoother regions (away from the centre in the Aig case) the shapes at the 
two /5's should be more similar to each other. This is indeed observed; in the Eu 
case, the largest differences are at Rt ~ 2 or 3 instead of Rt = as for the Aig 
symmetry, because the distributions peak at these values. This effect has not 
been mentioned in earlier works ^. 

At flrst sight, it appears that the only cases where the results for the two /9's are 
consistent are S{Aig) and EiiAig), whereas in the 2-d sums of Figs. 4-6 other 
cases, such as S{Eu), El{Eu) and even Si^A'^g), ET{A[g) seem to show reasonable 
scaling. The possible reasons for this are twofold: 

i) Figs. 7-9 only show the proflle in a single direction - along a lattice axis - 
whereas Figs. 4-6 are an average over all directions in a plane. In particular, this 
could have an effect on small values of Rt, where rotational invariance is most 
violated. 

ii) The curves depicted in Figs. 7-9 must be multiplied by a phase space factor 
2nRTdRT, when their contributions to any 2-d sum rule are estimated. Therefore, 
the values near Rt ~ get drastically reduced and, in the Aig and A[g cases, it 
is this region of Rt that is varying the most with (3. 

Assuming some function describing the continuum density we could apply a dis- 
cretisation procedure, e.g. simply averaging over cubes of volume a^, that simu- 
lates the flattening of the peaks in our flnite-a simulations. The latticized con- 
tinuum function could then be fltted to measured points. When this procedure 
is applied at both simulated values of the lattice spacing, we would get two cor- 
rected parameterisations. If these two agree this then would suggest that this 
transverse distribution would apply in the continuum limit. This would be the 
way to compare with continuum models of the flux tube: indeed we flnd qual- 
itative agreement with the dual QCD model of Refs. |Ty, |TT] for the Aig 
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profiles. 

An interesting feature in the A[g profiles is a local minimum (a dip) outside the 
centre of the tube as predicted by the N=l Isgur-Paton model ||12| for the energy 
density. In Fig. ^ the action can be seen to have a plateau at 3 < Rt < 5 unlike 
in the Aig case. For El no evidence of a dip is found, whereas for Et the data 
with f3 = 2.4 hints at a minimum for Rt ~ 3 with the 4 < Rt < 7 values being 
above zero unlike the Aig case, which is consistent with zero for Rt > 4. 

A better statistical accuracy is achieved in the transverse profiles of fiux-tubes 
with interquark separation i? = 4, 6 at /3 = 2.4, 2.5 respectively, corresponding 
to interquark distances of 0.473(5) and 0.503(7) fm. These are shown in Fig. |10| 
with the longitudinal and transverse components of the action plotted separately. 
There is now a clear dip in the transverse action profile, again at Rt ~ 3, with 
a corresponding maximum at Rt ~ 5. A similar dip and maximum are seen 
in the longitudinal action at (3 = 2.5, whereas the f3 = 2.4 data has a plateau 
at 3 < Rt < 6. As in Fig. |^, no evidence for a dip is found for E^, whereas 
the Et data shows a clear dip at /3 = 2.5 with a minimum at Rt = 4. At 
f3 = 2.4 there is no minimum but the decay as a function of Rt is slow, the value 
at Rt ~ 6 being an order of magnitude higher than in the corresponding Aig 
case. However, again it should be emphasized that Table 2 indicates for the A[g 
state considerable contamination from neighbouring states. Therefore, any nodal 
structure possibly present in a pure A[g state could well be smoothed out by 
interference effects. Also, it should be added, that the Isgur-Paton model which 
suggests such dips is less applicable for these smaller values of R. 



4 Determination of /3-functions from sum rules 



In Ref. 0, by imposing the condition on the interquark potential V that 
dV{R)/da\Fi = 0, the following three sum rules were derived relating V to spatial 
sums of the electric and magnetic colour. 

~^ ^y + R%]+So = J2S = -T.i^L + 2ST + 2BT + BL) (14) 



b V OR 

1 



^(^ + ^^)+^o = J:El = U-^l + Bl) (15) 

= T.ET = T.i-^T + BT). (16) 



Here the generalised /^-functions are defined considering an asymmetric lattice. In 
the notation of Ref. §, they are b = d/3/d\na = 2{S + U) and f = {U-S)/{2(3). 
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Figure 7: The colour flux contributions corresponding to the action [S), longi- 
tudinal (El) and transverse energy (Et) sum rules of Eqs. [l^-[l6| for the static 
quark potential. These are shown in units of a(2.5) versus transverse distance 
Rt at the mid-point {Rl = -R/2) for separation i? = 8, 12 at /3 = 2.4, 2.5. The 
data are for the symmetric ground state (Ai^ representation). 



In Eqs. p!4|-p!q, 5*0 and Eq are the self-action and -energy associated with the 
quarks and are, therefore, independent of R. The same self-energy Eq is expected 
for both orientations [L and T) of the colour electric field. 

The three-loop perturbative expression for b in terms of a = g'^/^Air) = l/(7r/3) is 
b = -0.37151(1 + 0.49193a - 0.9795^2 + ...). On the other hand, for / we have 
/ = 1 — 0.4567ra; — 0.257rQ;6 + . . . where we can insert the expression for b giving 
/= 1 - 1. 1408a + .... 

The aim is to now extract estimates of b and / in the non-perturbative situation 
encountered in practice. This will be carried out in various ways, each of which 
has its own advantages and disadvantages. However, their common feature is 
that, on the LHS of each sum rule, the potential V is measured on the lattice 
using Wilson loops W{R,T) and, on the RHS, as discussed in the introduction. 
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o o |3=2.4 Eu 

• • |3=2.5 E„ 




Figure 8: As Fig. |^ but for the first gluonic excitation {E^ representation). 



the colour fields S and S at a point r, are measured using Eq. |I| involving the 
same loops. Unfortunately, this strategy is complicated by two features in Eqs. 
n-TB. Firstly, the self-energies are unknown and so their effect must be removed 
by considering differences between the equations for different values of i? - or for 
different gluonic states Aig, E^, ■ ■ ■ with the same R, a possibility not considered 
here. Secondly, each equation contains dV/dR. Even though, the potential V 
itself can be readily determined as a by-product of the B measurement on a 
lattice, to determine dV/ OR introduces some uncertainty. All of the estimates in 
this section are made using the Aig ground state, where we have the best signal. 
The other gluonic states have such large errors that sensible values of 6, / can 
not be extracted. 



4.1 Method 1: Fitting the sum rules 

The most direct approach is to measure £ and B over all space and to then 
perform the spatial sum giving the RHS of Eqs. ^^-|T^. In practice, "spatial 
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o |3=2.4 A',„ 
• |3=2.5 A',° 




Figure 9: As Fig. |^ but for the second gluonic excitation {A[g representation). 



sum" means a sum over a lattice that has a hnear size twice that of the maximum 
R considered i.e. for (3 = 2.4 (2.5) up to i? = 8 (12) on lattices with spatial size 
16^ (24^). With V known numerically from Eq. the derivative can also be 
estimated. The function b is then obtained using Eq. |T^ by plotting J^i^L + 
2£t + Bl + 2Bt) vs. V{R) + i?^^^ and performing a hnear fit, as shown in 
Fig. |TT[ The R = 1 points were not included in the fits due to the artefacts they 
contain, while the i? = 12 point at /3 = 2.5 was excluded because of its significant 
excited state contamination (see Table There are four sets of data for each 
l3 corresponding to the correlation of the sum over electric and magnetic fields 
taken at time intervals T = 3, . . . , 6. The (3 = 2.4 data was scaled by multiplying 
with 2.4/(2.5p) to have the same units as the /5 = 2.5 plots. The results of the 
fit in Fig. |ll] can be read from the second column of Table ^. 

For (3 = 2.4 the function b has reached a plateau at T = 5, giving a best estimate 
of —0.312(15), whereas for j3 = 2.5 a plateau has been reached only at T = 6 
with —0.323(9) being our best estimate. The self-energy estimates 5*0 are also 
seen to reach plateau at —1.2(1) and —1.5(1) for (3 = 2.4, 2.5 respectively. 
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Figure 10: As Fig. |^ but with longitudinal and transverse components of the 
action presented separately, and for interquark separation R = 4,6 at P = 2.4, 2.5 
respectively. 



Even though the P = 2.4 and 2.5 lines in Fig. |TT] are almost parallel as in Fig. |l], 
they are separated by the difference in their self-energies. The strategy of fitting 
the data with a straight line is effectively taking differences of Eq. [1^ evaluated 
for different values of R and so avoids the need to know 5*0 explicitly. 

In principle, the function / can be extracted from either of the two sum rules in 
Eqs. |15| and ^ However, as shown in the fourth and sixth columns of Table ^ 
these predictions exhibit much more variation and have much greater errors than 
those for b. The main reason for this is due to the larger variation of the sums 
of the differences of electric and magnetic fields. In particular, the values of /(/) 
at /5 = 2.4, T = 5, 6 and (3 = 2.5, T = 4, 5, 6 are essentially undetermined. 

The self-energies Eq are not consistent with zero for /(/) at /? = 2.5 as they are 
for /(J, II) at /5 = 2.4 and /(//) &t (3 = 2.5. As Eq is the difference of: 

1) Vo/(4/5/) coming from the two-body potentials in Eq. |^ and 

2) the self-energy in E^, Et for /(J, II) respectively. 
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Figure 11: Data corresponding to Eq. with best linear fits. 

we can see that these two cancel for all cases except /(/) at /5 = 2.5. As the 
self-energies should be isotropic, this is probably caused by inaccuracies in deter- 
mining 

This inaccuracy in determining / can be greatly reduced if the sum rules in eqs. 



15 and 16 are fitted together. At /? = 2.4, 2.5 data for the former sum rule is 



taken at T = 3, 4; T = 3 respectively, while for the latter it can be taken at any T 



from 3 to 6. The function / obtained in this manner is presented in Fig. and 
Table ^ and can be seen to lead to a much more accurate estimate of /. Our best 
estimates are 0.65(1) and 0.68(1) at = 2.4, 2.5 respectively. Also the values of 
self-energy are now more stable at Eq 0.01(1) - a number that is about two 
orders of magnitude smaller than the self-actions 5*0- 

In the above, the parameter p has been extracted and found to be 1.411(13). 
This then suggests as a direct estimate of b, averaged over the f3 range of 2.4 to 
2.5, the value 

Ad 2.5-2.4 
^ -0.290(8). 



Aln[a] ln[a(2.5)] - ln[a(2.4)] 
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Table 3: Estimates for b and f at f3 = 2.4, 2.5 at different T values. 





T 


b 


5*0 


/(I) 


Eq 


/(II) 


Eq 


2.4 


3 


-0.357(3) 


-1.04(2) 


0.63(3) 


0.004(4) 


0.63(5) 


0.005(3) 




4 


-0.336(5) 


-1.11(3) 


0.71(6) 


0.010(7) 


0.74(10) 


0.007(4) 




5 


-0.312(15) 


-1.24(10) 


1.5(1.3) 


0.06(2) 


0.84(28) 


0.01(1) 




6 


-0.317(21) 


-1.21(14) 


4.5(7.0) 


0.09(4) 


1.4(8) 


0.02(2) 


2.5 


3 


-0.389(4) 


-1.17(2) 


0.89(7) 


0.033(7) 


0.59(3) 


0.002(3) 




4 


-0.354(5) 


-1.32(3) 


1.41(23) 


0.059(10) 


0.64(5) 


0.003(5) 




5 


-0.333(7) 


-1.44(4) 


1.49(38) 


0.063(16) 


0.75(11) 


0.011(8) 




6 


-0.323(9) 


-1.49(7) 


6(11) 


0.10(3) 


0.74(16) 


0.009(12) 



Table 4: Combined fits of / at /3 = 2.4, 2.5 at different T values. The first T 
value refers to the data used for Eq. |T5|, the other to the data used for Eq. |T6[ 





T 




Eo 


2.4 


3,3 


0.647(7) 


0.006(1) 




3,4 


0.616(8) 


0.0008(10) 




3,5 


0.627(17) 


0.003(3) 




3,6 


0.612(18) 


0.013(4) 




4,3 


0.705(13) 


0.009(1) 




4,4 


0.671(13) 


0.004(1) 




4,5 


0.684(23) 


0.040(5) 




4,6 


0.659(27) 


0.030(6) 


2.5 


3,3 


0.694(11) 


0.010(1) 




3,4 


0.667(11) 


0.005(1) 




3,5 


0.682(16) 


0.008(2) 




3,6 


0.688(21) 


0.009(3) 
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Figure 12: Determination of f{I + II) by combining two sum rules. 



Even though this is admittedly a very crude estimate, it should represent the 
average value of h over this range of 13. However, it appears to be slightly smaller 
in magnitude than the average of the above estimates —0.312(15) at /9(2.4) and 
—0.323(9) at /9(2.5). The origin of this two sigma difference is not clear. 



4.2 Method 2: Combining the sum rules 

As shown in Ref. 0], one way to avoid estimating dV{R)/dR and S'o, Eq is to 
explicitly eliminate them by writing down Eqs. [T^[T6| for two different values of 
R. In this way 

, ^ 2(nfi.)-'^(fi.))(i + gS;:ggl!g )"' ^^^^ 

^ 2(3[nET)R,-nET)R, + nEL)R,-nEL)R,\ ^ ^ 
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At first sight this appears to be what is needed - expressions that involve quanti- 
ties that can be measured directly. However, in practice, there is a problem - b be- 
comes dependent on the differences J2{El)r-^ — J2{El)r^ and J2{Et)rj^ — J2{Et)r2 
from Eqs. p!5Hl6|, and, as seen above, the values of these differences are less ac- 
curate than J2 Sr^ — J2 Sr^ . The outcome of this strategy is given in Fig. |13 



There it is seen that the 6(2.4) results are consistent with those given by method 
1 - but have much larger error bars. However, the 6(2.5) results are essentially 
inconsistent with method 1. A similar problem arises with the values of / from 
Eq. Again /(2.4) is consistent with the earlier estimates of method 1 in Ref. 
[0 - but with much larger error bars. For example, with Ri,R2 = 2,6 we get 
/(T = 4,5) = 0.64(7), 0.64(15). However, compared with method 1, /(2.5) is 
badly out - rising to ~ 0.9(2) at T ^ 4, 5. It should be added that this is not a 
problem of the measurements being poorly distributed, since plotting the boot- 
strap values of b shows that the errors are not underestimates due to asymmetric 
non-gaussian bootstrap distributions. 
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Figure 13: Estimates of b from Eq. ^ for a) /? = 2.4 h) P = 2.5. 
This difference between the two estimates of b for (3 = 2.5 but not 2.4 is caused by 
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the inaccurate determination of the longitudinal energy sums I](-^l)-Ri 
in Eq. |l^ As can be seen in Fig. 0, for the J^El curves at /? = 2.5 the slope 
is rather erratic and indicates a much smaller value than at (3 = 2.4. This is 
also reflected in Table ^, where the /(/) value obtained using the J2 El values is 
already at T = 3 unrealistically high for (3 = 2.5, getting worse with increasing 
T. This means that the Z](-El)ri — Z](-Sl)/?2 underestimated. For J2Et the 
situation is more consistent, which can be seen in the larger slope in Fig. ^ and the 
reasonable behavior for the /(//) in Table | at /3 = 2.5. Thus J2{El) ji-^—J^iEi) 
is too small and J2{Et)ri — S(-E't)_R2 realistic, leading to underestimates of 6(2.5) 
using Eq. and an overestimate of /(2.5) from Eq. 0. The signal being worse 
for El than Et is somewhat surprising, since it is Et where the self-energy 
completely dominates. Therefore, one would have expected it to be harder to get 
a signal for Et, because it requires a more delicate cancellation of the self-energy. 

Why doesn't this happen at /? = 2.4? From Figs. |^ and ^ we can see that the 
J^El slope is larger and the J2Et slope slightly smaller than at j3 = 2.5. This 
is again reflected in Table where /(J) again gets larger with increasing T, but 
not as much as for (3 = 2.5 and with larger errors making the estimates consistent 
with a realistic value. At the same time, unlike for f3 = 2.5, f{II) also increases. 
This fortuitously leads to a realistic ratio of the differences of the sums in Eq. |1^ 
and a value of h consistent with method 1. 



4.3 Method 3: The large R Umit 

In the above, the derivatives of V are calculated numerically from the lattice 
form of the interquark potential in Eq. |^. However, for sufficiently large R {R> 
2), the continuum form of V (i.e. with [1/-R]l replaced with 1/-R) is a good 
approximation. When, in addition to this, the effect of the self-energies is removed 
by evaluating the sum rules at two values of R, Eqs. p!^p!6| reduce to 

. _ -2bs{R^-R,) ^^g^ 



J2 



R2 



^ 2f3[E{EL)R.-E{EL)R,] ^ ' 

fill) = - .21) 

^ 2p[nET)R,-nET)R,\ ^ ^ 



As seen in figure O a) for (3 = 2.5, Eq. O gives estimates in agreement with the 



fits in Table ^ for Ri, R2 = 2, . . . ,12. However, this is not surprising, since there 
the results are an average over a range of R values, whereas Eqs. |T^-^ can be 



considered as an average using simply two values of R. A similar situation holds 
for (3 = 2.4. 
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Figure 14: Estimates of b at j3 = 2.5 using a) the continuum form of V{R) 
(Eq. |19|) and, additionally, b) a constant longitudinal profile approximation for 
the flux tube with e = 0.5 (Eq. g2[) 

Of more interest is the large R limit, where it is expected from string models 
that both the action and energy fiux-tubes should have a form that is essentially 
constant for all Rl in the range —R/2 < Rl < R/2 with i?-independent self- 
action or -energy terms concentrated at Rl ~ ±R/2. As can be seen in Fig. I, the 
sums over the transverse plane in the middle of the quarks agree at /3 = 2.4 and 
2.5, supporting the accuracy of our estimates in this region. Therefore, the action 
difference in Eq. |I9], AJ^S = (J2 Sr^ — J2 Sr^) should be well approximated by 
A E ^ = [RiSr, {Rl = 0)-i?2^i?2 (Rl = 0)] and similarly for A E and A E 
in Eqs. |2g and ^. 

In practice, since Ri^2 are not very large, it is probably more realistic to include 
a correction for the 'overshoot' of the plateau term at Rl ~ ±i?/2 i.e. 

AJ2S^AJ2S'= [R,{eS)SR,{RL = 0) - R2{eS)SR,{RL = 0)], (22) 

where Ri{eS) = Ri + 2e. Here e is the 'overshoot' at each end of the plateau. 
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From Figs, ^j-^, it is seen that e must be approximately one half of a lattice 
spacing. 

The results at (3 = 2.5 with e = 0.5 are shown in Fig. |T^ b). The values of 
b approach those given by method 1 as R12 increases, agreeing when R12 = 
6, 12. The disagreement at the smaller i?'s is mainly due to significant self-energy 
contributions at the midpoint, which don't cancel as they vary with R. The good 
results from R = 6 thus show that the self-energy contributions at midpoint are 
negligible at these interquark separations. The same approximation can be used 
in the spatial sums of method 1. With R = 2,3 excluded from the fit of Eq. |14| 
we get b = -0.35(4) at T = 6 and /5 = 2.5. 

The extracted value of b is only weakly dependent on e even for the rather small 
values of Ri^2 used here. This situation will only improve as Ri^2 increases. In 
principle, this method has two advantages over the earlier ways of extracting b 
and /: 

1) It avoids the need for any delicate numerical cancellations of self-energies, since 
they are assumed to have canceled exactly in the A's. 

2) Only the single 2-dimensional sum over the Rl = plane is necessary. This 
avoids the full 3-dimensional sum of Eqs. |l^ ^T^ . This estimate can be further 
improved by averaging the 2-dimensional integrals over the smallest values of 
Rl. However, in practice this would not lead to savings in computer time, since 
averaging over all positions and orientations of the Wilson loops would, in any 
case, mean using all possible planes in the lattice. 

It has one disadvantage, however, in that it assumes a string-like longitudinal 
dependence of the colour flux distribution. 

For large R one would expect that the string tension is given by the longitudinal 
energy density in the transverse plane at the midpoint: 

bs = 2(3fEL{RL = 0). 

Taking the values of El{Rl = 0) from Fig. | and using our best estimates of / 
gives ^bs{2A) = 483(15), ^6^(2.5) = 477(24) MeV. These are close to the values 
given by fits to experimental spectra. 

4.4 Comparison with other approaches 

There are two main ways to extract /5-functions in lattice gauge theories. First, 
one can measure observables at different values of some parameter (such as a cou- 
pling or a quark mass) and then try to estimate the response of the observable to 
a change in that parameter. This can be carried out using either finite differences 
or an interpolating function. Both of these introduce a systematic error either 
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from the use of a finite interval or from the choice of the interpolating function. 
In the second method, one can use the fact that response functions are related 
to correlation functions via lattice sum rules. This method has less systematic 
errors but often leads to large statistical errors, because the observables involve 
delicate cancellations. 

In the case of SU(2) gauge theory, the first method has been used recently in Ref. 
[0, where b was estimated by measuring the critical temperature at six values 
of the coupling and using an exponential ansatz or a spline interpolation for the 
P dependence of a. In the same work, non-perturbative derivatives of couplings 
in the time and space directions with respect to the asymmetry factor ^ = a^/at 
were estimated from h and measurement of the free energy. In Ref. lattice 
spacings were determined from the string tension at seven couplings in the range 
2.3 < (3 < 2.85 and fitted with an ansatz of a linear a dependence of Aiat- The 
errors on these estimates do not include systematic effects from the use of a linear 
approximation. In Ref. the lattice spacing was determined at six couplings 
ranging from 2.3 < (3 < 2.55 through Sommer's equation rQF(ro) = c discussed 
after Eq. ^. Fitting of the three-loop relationship between a and /3, with extra 
terms and systematic errors estimated using different values of c and forms of the 
fitted function, lead to results consistent with Ref. iSll. 



The sum rule method was first explored in Ref. ||T3|, where dV{R)/d(3 was fitted 
from the time dependence of Eq. [I| and b then extracted through a second fit of 
dV{R)/dj3 vs R. Optimistically also potentials at i? = 1 were included in the fits. 
As this was only a feasibility study the estimate of 6 = —0.25(2) at /5 = 2.5 should 
be considered as preliminary. A direct formula for b, Eq. involving differences 
of two-body potentials at two values of R and their action sums was derived in 
Ref. An average over pairs of R values gave the estimate b = —0.35(2) at 
(3 = 2.4. The values at some R pairs were outside this estimate and doubts 
about the convergence in T remained, as the correlations were calculated only up 
to T = 4 for larger i?'s. 

Table ^ collects the most reliable estimates of 6 at /3 = 2.4 and 2.5, whereas 
the available estimates of / are shown in Table |^. Here it should be noted 
that an effective coupling larger than the bare coupling moves the perturbative 
estimates of / closer to the non-perturbative values, whereas the opposite happens 
for b. A perturbative evaluation of b is thus unreliable. The non-perturbative 
estimates of b agree with each other, except for the one from Ref. 0, which 
has a significant systematic error. Averaging the three consistent estimates gives 
6(2.4) = -0.306(6), 6(2.5) = -0.316(4), which are 78%, 81%, respectively, of the 
three-loop predictions. For / our estimates are not far from those in Ref. 
with similar ratios of the non-perturbative and perturbative estimates as for b. 
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Our estimate 


Ref. i 


Ref. i 


Ref. i 


3-loop PT 


2.4 
2.5 


-0.312(15) 
-0.323(9) 


-0.3018 
-0.3115 


-0.305(6) 
-0.312(2) 


-0.330(4) 
-0.340(4) 


-0.3893 
-0.3889 



Table 5: Comparison between values of 6 = (9/5/9 In a 



p 


Our estimate 


Ref. i 


PT 


2.4 


0.65(1) 


0.72 


0.85 


2.5 


0.68(1) 


0.74 


0.86 



Table 6: Comparison between values of / = ([/ — S)/{2(3) 



5 Conclusions 

In this paper the spatial distribution and nature of the electric and magnetic 
colour fields between two quarks are measured for (3 = 2.4 and 2.5. We discuss 
carefully the theoretical expectations for scaling versus (3 of these distributions 
and compare with our results. For the observables with a well-controlled contin- 
uum limit (three-dimensional sums), scaling is investigated and found to be good 
in most cases of interest. For more differential observables (transverse sums and 
-profiles) the changes seen between the two values of P can be explained quali- 
tatively from the discretisation: the plateaux stay the same, while the peaks get 
higher for smaller a. Self-energies, as measured by transverse plane sums, diverge 
in a manner suggested by leading order perturbation theory, with the divergence 
being faster in the El case because of the longitudinal extent of the plaquettes 
measuring the dominant electric field. By comparing various combinations of the 
spatial sums of the £, B fields with the interquark potential V{R) ± dV{R)/dR, 
estimates can be made of the generalised /3-functions b and / - see Eqs. p!4Hl^. 

There are two problems that prevent the direct use of these equations. Firstly, 
since the whole calculation is performed on a lattice, V{R) is only known at 
discrete values of R, so that the values of dV{R)/dR can not be measured di- 
rectly. Secondly, both V{R) and the spatial sums have self-action or -energy 
contributions. To minimize or avoid these problems the extraction of b and / 
can be carried out in different ways each of which has its own advantages and 
disadvantages concerning statistical and systematical errors. 

The most direct approach - method 1 of subsection |4. Ij - is to parameterize V{R) 
so that dV{R)/dR can be readily calculated. This overcomes the dV{R)/dR 
problem at the expense of some systematic error introduced by the form of pa- 
rameterisation of V{R) used. By plotting, as a function of R, V{R) ±dV{R)/dR 



28 



versus the various spatial sums, the slopes give immediately b and /. This re- 
sulted in 6(2.4) = -0.312(15) and 6(2.5) = -0.323(9). However, the value of / 
was much more uncertain when the two sum rules - one involving J2 El and the 
other J2 Et - were used separately. When these two sum rules were fitted simulta- 
neously, a better result emerged, /(2.4) = 0.65(1) and /(2.5) = 0.68(1). Already 
at this stage, it could be seen that the data causing most of these uncertainties 
are those involving the longitudinal energy J^El especially for (3 = 2.5. 

One way to avoid the dV{R)/dR and self-action and -energy problems is to 
combine the original three sum rules in such a way as to eliminate the problems 
- see Eqs. p!7|-p!^ for method 2. However, in spite of the ideal form of these 
equations they result in poor estimates of both b and /. For /3 = 2.4, b and / 
agree with method 1 but with large statistical errors - see Fig. |T3| a). But for 
(3 = 2.5 the values of b tend to be too low and of / too high - both with large 
error bars. These differences can be directly attributed to the inaccuracy of the 
data for J2 El- 

The third method for extracting b and / exploits the form expected of V{R) 
and the flux tube proflles in the limit R ^ oo from string models. When a 
continuum parameterisation of V{R) is used and two R values are subtracted, 
the sum rules reduce to those in Eqs. p!^-|2T| and these yield results consistent 
with those of method 1, but with larger statistical errors. In the large R limit, a 
further approximation is to replace the flux tube proflle by one that is constant 
between the quarks (— -| < Rl < -f ) and dropping to zero rapidly for Rl > ^ 
and Rl < ""f • At the largest i?'s the results are consistent with method 1, but 
with large statistical errors. These results are only very weakly dependent on the 
point where the proflle drops to zero beyond the positions of the quarks. 

As seen in tables ^ and our best estimates for b and / agree with other recent 
non-perturbative estimates, most importantly the flnite temperature approach of 
Ref. 0. It thus seems safe to conclude that order effects in the extraction of the 
/3-function are small at the /5-values studied using the methods described. Thus 
we have a unique non-perturbative /3-function which describes the deviations from 
asymptotic scaling at these values of the coupling. 

For the state with A[g symmetry, our data shows the existence of a dip in the 
action and transverse energy distributions away from the centre of the flux-tube 
in the transverse plane between the quarks. This qualitatively conflrms the pre- 
diction of the Isgur-Paton flux-tube model for the energy distribution. No such 
node is seen for the longitudinal energy, and for the action it is more pronounced 
in the transverse component. 

In the future we plan to study four-body flux distributions and their relationship 
to two-body flux-tubes. Sum rules will be used verify the correctness of the 
measured distributions. 
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